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Based on first-principles calculations, we predict that the magnetic anisotropy energy (MAE) 
of Co-doped Ti02 sensitively depends on carrier accumulation. This magnetoelectric phenomenon 
provides a promising route to directly manipulate the magnetization direction of diluted magnetic 
semiconductor by external electric- fields. We calculate the band structures and reveal the origin of 
carrier-dependent MAE in /c-space. In fact, the carrier accumulation shifts the Fermi energy and 
regulates the competing contributions to MAE. The first-principles calculations provide a straight- 
forward way to design spintronics materials with electrically controllable spin direction. 
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Diluted magnetic semiconductor (DMS), exhibiting 
both ferromagnetism and semiconducting properties, has 
been considered to be a promising candidate of spintron- 
ics. It provides a potential route to add the spin de- 
gree of freedom to conventional charge-based electronic 
devices, for example, adding magnetic recording capa- 
bility to current semiconductor information processing 
unit. However, there are two major challenges, one is 
how to raise the Curie temperature (Xc), and the other 
is how to control the magnetization direction. Since the 
electric-field-induced room-temperature ferromagnetism 
has been demonstrated in cobalt-doped titanium dioxide 
(CoiTiC^)^ where the ferromagnetic exchange is consid- 
ered to be mediated by carriers and tuned by electric- 
field, the left issue is how to manipulate the magnetiza- 
tion direction. 

Conventionally, the magnetization direction is con- 
trolled by magnetic-field, which is unsuitable for 
ultrahigh-density magnetic storage and integration with 
electronic devices. Therefore, explorations aiming at ma- 
nipulation of magnetization direction directly by electric- 
field have emerged. Early experiment 2 reveals that the 
coercive force Hq depends on the variation of carrier den- 
sity managed by electric-field in Mn-doped InAs. This 
phenomenon hints that there is a correlation between the 
magnetic anisotropy and carrier density in DMS. Further 
work 3 exhibits the rotation of magnetization direction by 
applying electric-field on Mn-doped GaAs, showing the 
direct connection between magnetic anisotropy and car- 
rier density. 

When an external electric-field is applied, there will be 
an accumulation of electrons or holes, yielding a shift of 
Fermi energy. Since magnetic anisotropy energy (MAE) 
is mainly determined by the band structure near the 
Fermi energy, the shift will impact the MAE. Early first- 
principles calculations of MAE in transition metal bulks 
and thin films have shown this effect by the so-called 
electron-filling technique, where electrons are added to 
or removed from the system under studyj 4 -"— We expect 



that carrier accumulation may also impact the MAE in 
DMS. Moreover, it will be easier to realize a carrier accu- 
mulation in DMS in experiment than magnetic metals. 

In this letter, for Co:Ti02, a typical DMS, our numer- 
ical calculations demonstrate that the carrier accumula- 
tion yields an oscillation of MAE. Instead of analyzing 
MAE within traditional single-ion anisotropy theory, a 
straightforward analysis based on detailed band struc- 
ture near the Fermi energy is proposed to explain the 
dependence of MAE on carrier accumulation. 

Based on density functional theory (DFT), first- 
principles calculations were carried out on Co-doped 
Ti02 anatase using Perdew-Burke-Ernzerhof (PBE) 
parameterization^ of generalized gradient approximation 
(GGA) as implemented in VASP package.- The primi- 
tive anatase cell was fully optimized. Then, a 2 x 2 x 1 
supercell was created with one Ti atom substitued by Co 
(Coo.o625Tio.9375 02), and the atomic positions were al- 
lowed to relax. After the optimization, GGA+ U (GGA 
plus on-site Coulomb repulsion) approach 9 was employed 
in electronic structure calculations. We applied extra 
Coulomb repulsion to Ti-d orbital (2 eV)^ and Co-d or- 
bital (2 eV)il. The plane wave cut-off energy was 500 eV. 
The tetrahedron method with a 5 x 5 x 4 A;-mesh grid was 
employed for the integration in Brillouin zone. The ac- 
curacy of electronic iterations was up to 10 -6 eV. The 
MAE, merely considering the contribution from spin- 
orbit coupling, was calculated as MAE = i£[ioo] — -^[ooi]? 
where 22[ioo] an d ^[ooi] were the total energy with mag- 
netization directions along [100] and [001], respectively.— 
Carrier accumulation was simulated by modifying the to- 
tal number of electrons, assuming a homogeneous back- 
ground charge. The MAE is calculated for multiple 
charge configurations. 

FIG. Q] shows the MAE and the magnetic moment as a 
function of the electron number added (SN), where neg- 
ative 5N means that electrons are removed. It is obvious 
that both the MAE and the magnetic moment sensitively 
depend on electron- filling. In neutral Co:Ti02, the MAE 
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FIG. 1. (Color online) (a) The total magnetic moment, mag- 
netic moment on Co atom and (b) MAE of Coo. 0625^0.937502 
depend on electron-filling. 



is negative, implying the anatase c-plane is the easy plane 
of the magnetization. However, when electrons are re- 
moved, the easy-axis rotates out of the c-plane, and the 
c-axis becomes the easy-axis, and the magnitude of MAE 
increases steeply. In addition, when electrons are added, 
the magnitude of easy plane magnetic anisotropy begins 
to reduce and then becomes zero at SN = 1.0 and 1.5. 
With more electrons doped, the easy-axis aligns along 
the c-axis, but the magnitude of MAE is weak, less than 
0.1 meV. According to these results, the magnetization 
direction will switch from in-plane to perpendicular by 
varying carrier density, which can be achieved by carrier 
accumulation and tuned by external electric-field. 

When SN decreases from to -2.0, the total magnetic 
moment of the super cell (M to t) linearly increases from 
1.0 to 3.0 /iB, however, the local magnetic moment on 
Co (Mco) remains steady. In contrast, when SN varies 
from to 1.0, both M tot and Mq decrease to /iB, in- 
dicating the system becomes non-magnetic. The nearly 
coincidence of the two curves suggests the magnetic mo- 
ment mainly contributes from the Co atom at this range 
of SN. Further, we continue to increase SN from 1.5 
to 3.0, the system regains magnetic moments. Finally, 
when more than three electrons per supercell are added, 
the M to t drastically rises to 2.0 /iB and then linearly in- 
creases to 2.5 /iB, while Mq levels off. 

To understand the impact of electron-filling on mag- 
netic moment and MAE, we first check the electronic 
structure of Co:Ti0 2 . The density of states (FIG. [2j) for 
neutral Co:Ti02 with GGA+ U shows that the band gap 
of pristine Ti02 is preserved, with the impurity bands 
lying in the gap. The well isolated impurity bands are 
mainly from the Co d-orbitals, which hybrid with O p- 
orbitals near the valence band maximum (VBM), and 
lead to slight magnetization of O atom as shown in the 
spin density map (FIG. [2j). The majority of Co £2 man- 
ifold is completely occupied, while the minority of £2 
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FIG. 2. (Color online) Spin density maps (inset) and density 
of states for Coo.o625Tio.937s02 with GGA+ U (Uri-d = 2 eV, 
Uco-d = 2 eV) approach, the vertical dash line at eV is the 
Fermi energy. Solid line plot, total; red filled plot, Co-d £2; 
blue filled plot, Co-d e, respectively. 



manifold is split into an occupied doublet (d xz?yz ) an d 
an empty singlet (d xy ). The insulator ground state is 
different from the half-metallic ground state predicted 
by GGA*i2r— The total magnetic moment of the cell 
(M t ot = 1-0 ^b) and the spin density map (FIG. [2]) of Co 
atom with a shape of d xy orbital suggest a low-spin state, 
(^l) 3 ^) 2 ' which is consistent with previous beyond-DFT 
results J^ii^ 

The electron- filling causes an obvious shift of the Fermi 
energy (FIG. [3]), if the projected density of states with 
different SN is aligned to the deep O s-orbital. When 
one electron is added (SN = 1.0), the singlet (d X y) i n 
minority spin is occupied, and the system becomes non- 
magnetic. In addition, the Co-d e manifold is pushed 
into the conduction band (CB) of Ti02 host. When two 
electrons are added, the Fermi energy shifts further into 
the host CB, and results in the so-called "band- filling ef- 
fect", i.e., the host conduction band minimum (CBM) 
will be first occupied and then the e manifold of Co 
will be partially filled. 16 Similarly, when electrons are 
removed, they will be partially removed from the Co t^. 
Therefore, Mq almost maintains a constant value with 
SN < and SN > 1.0 [FIG. da)], which is called the 
negative-feedback charge regulation in Ref. [I?]- 

It is well-known that the spin-orbit coupling (SOC) in- 
teraction of 3d transition metal elements is much weaker 
than the crystal- field split, and that MAE can be esti- 
mated by single- ion anisotropy theory j 18 i 19 However, the 
system becomes metallic when SN < and SN > 1.0, 
and consequently, the charge state of Co is no longer 
well-defined under electron-filling. Then, the traditional 
single-ion anisotropy theory, where a well-defined charge 
state of magnetic ion is required, might not be suitable 
under this circumstance. 
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Projected Density of States (states/eV) 

FIG. 3. (Color online) The projected density of states 
(PDOS) of Coo.o625Ti .93750 2 with different SN. The PDOS 
is aligned by the deep 0-5 ( ~ —17 eV). Therefore, the 
electron-filling causes an obvious shift of the Fermi energy 
(horizontal wine line). Gray rilled plot, total; magenta rilled 
plot, Ti-d; blue filled plot, Co-d e; red filled plot, Co-d ti. 
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respectively. It is obvious that the manifold splits into 
a singlet (d xy ) and a doublet (d xz , yz ), due to the local 
D2d symmetry of the Co dopant. The doublet further 
splits along some directions, e.g. Z to R, due to the dis- 
persion, when the translation symmetry is considered. 
As a result, the contribution to MAE can be divided 
into three categories: (i) the SOC interaction between 
occupied <i xz , yz and empty d xy , (ii) the SOC interaction 
between occupied d xz (d yz ) and empty d yz (d xz ), and (hi) 
the SOC interaction of the degenerated states inside the 
doublet d xz ,yz- The sign of MAE can be estimated by 
summing above three contributions. 

For (i) the SOC interaction between occupied <i xz?yz 
and empty d xy that have different magnetic quantum 
numbers, the perturbation is through L x operator and 
yields a negative contribution. 6 - For (ii) the SOC interac- 
tion between occupied d xz (d yz ) and empty d yz (d xz ) that 
share the same magnetic quantum numbers, the per- 
turbation is through L z operator and yields a positive 
contribution. 6 * Note that the energy split between d xz and 
d yz is almost tenth of that between d xz , yz and d xy . The 
contribution from (ii) can be 10 times larger than (i). For 
(iii) degenerated state, L is unquenched. In the subspace 
spanned by \xz) and \yz), the SOC Hamiltonian can be 
written as 



U = \L-S 
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where is the angle between L and S. Our calculation 
shows that is also equal to the angle between spin di- 
rection and the c-axis, since L is always along the c-axis 
in the calculation. The eigenvalues of the Hamiltonian 



FIG. 4. (Color online) Band structure of minority spin near 
the Fermi energy (horizontal dash line at eV) along the edge 
of the irreducible Brillion zone (IBZ) without SOC, (a)SN = 
and (b) SN = -1.0. The Co-d bands are plotted in red. 



On the other hand, the MAE can be obtained by in- 
tegrating the net contributions of the SOC interaction 
between the 3d subbands in A;-space.— i 2Q i 21 The contri- 
butions from degenerated and non-degenerated pertur- 
bations results in the first-order and second-order con- 
tributions, respectively. For the non-degenerated part, 
the contribution to MAE depends on the interaction be- 
tween the occupied and empty states^ It should be noted 
that the degenerated contribution could be as important 
as the non-degenerated part, although the degeneracy 
occurs only in a small portion of the Brillouin zone^ 2 - 
Thus, we compare the unperturbed band structures near 
the Fermi energy for different SN. 

Due to the fully occupied ti manifold in majority spin 
and the low-spin configuration, the contribution from the 
spin- flip terms has been omitted ^ 2 1 i 23 We only plot band 
structures of minority spin near the Fermi energy for 
SN = and -1.0 in FIG. SJa) and (b) without SOC, 



£ 1 = B COS 0;£ o = -iM C o S 0. (2) 

If there is only one electron in these states, e.g. SN = 
— 1.0, the low- lying level Eo will be occupied. As a result, 
the energy of the system depends on with the minimum 
at = 0, i.e. the easy-axis is parallel to the c-axis. 
Therefore, we conclude that, for degenerated sates, the 
contribution to MAE is positive. 

With above analysis, the results of MAE for SN = 
and SN < can be explained qualitatively. For SN = 0, 
there is only perturbation of category (i), resulting in the 
moderate negative MAE. For SN = —1.0, there are all 
three contributions, where both categories (ii) and (iii) 
are positive contributions, and the magnitude of category 
(ii) is larger than (i). Thus, the MAE is positive. When 
more electrons are removed, the proportion of the posi- 
tive contribution from (ii) increases, and the magnitude 
of MAE is larger than that of SN = -1.0. 

The dependence of the total energy on the angle (6) 
between the quantum axis of spin and the anatase c-axis 
(FIG. [5} supports our argument quantitatively. The co- 
efficient of cos#, the first-order perturbation, indicates 
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FIG. 5. (Color online) The dependence of the total energy 
on the angle (0) between the quantum axis of spin and the 
anatase c-axis, (&)5N — and (b) SN = -2.0. 



the contribution from degenerated perturbation. For 
SN = 0, it is less than 10 -6 eV, implying the contri- 
bution from degenerated perturbation can be ignored. 
The sign of the coefficient of sin 2 the second order per- 
turbation, includes the competing contributions from the 
non-degenerated perturbation. 

For SN > 1.0, the MAE alters its sign to positive (the 
easy-axis parallels to the c-axis) and decreases its mag- 
nitude to one order smaller than others. As shown in 
FIG.[3j when SN > 1.0, the e manifold becomes partially 



occupied. In fact, there will be a strong Jahn- Teller dis- 
tortion for d 7 electronic configuration in low-spin state 
under octahedral crystal field. Consequently, the MAE 
will be weak, as Jahn- Teller effect increases the split be- 
tween occupied and empty states and lowers down the 
total energy in general. In fact, without constraining the 
atomic positions, the total energy is lowered by 1.028 eV. 

In summary, our first-principles calculations predict 
that the MAE in Co:Ti02 can be controlled by carrier 
accumulation. This magnetoelectric phenomenon in this 
DMS system allows an effective manipulation of the mag- 
netization direction directly by external electric-field or 
voltage. To interpret the behavior of MAE, the elec- 
tronic structures with electron-filling are examined. A 
self-regulated feedback effect of local magnetic moment 
on Co has been discovered. We propose an explanation 
for MAE for metallic state where the charge state is not 
well-defined. The MAE is discussed in &-space based on 
the band structure near the Fermi energy. Based on the 
perturbation method, the contribution to MAE has been 
divided into three categories. The shift of the Fermi en- 
ergy caused by electron- filling regulates the contributions 
from different categories, and consequently determines 
the sign and magnitude of MAE. 
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